Effects of a Cut, Lorentz-Boosted sky on the Angular Power Spectrum 
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The largest fluctuation in the observed CMB temperature field is the dipole, its origin 
being usually attributed to the Doppler Effect - the Earth's velocity with respect to the 
CMB rest frame. The lowest order boost correction to temperature multipolar coefficients 
appears only as a second order correction in the temperature power spectrum, Cg. Since 
v/c ~ f0~ 3 , this effect can be safely ignored when estimating cosmological parameters [1H7]. 
However, by cutting our galaxy from the CMB sky we induce large-angle anisotropics in the 
data. In this case, the corrections to the cut-sky Cis show up already at first order in the 
boost parameter. In this paper we investigate this issue and argue that this effect might 
turn out to be important when reconstructing the power spectrum from the cut-sky data. 

I. INTRODUCTION 

In the last decade, cosmology has become a precision science, driven largely by measurements 
of the fluctuations in the Cosmic Microwave Background (CMB) temperature. These fluctuations 
- nearly 10 5 times smaller than the average temperature of the universe - are the primary window 
onto most cosmological parameters [I]. Because the fluctuations are so small, one must know the 
contributions to the measured temperature field from foregrounds and other contaminants. A great 
deal of effort has been put in to characterizing foreground signals from dust, synchrotron, and free- 
free emission [2]. Even so, when inferring cosmological parameters significant fractions of the sky 
are usually omitted ("cut") from full-sky maps in order to minimize effects from non-primordial 
sources [3]. 

Surprisingly, one particular known systematic effect - the distortion of the CMB radiation due 
to our motion relative to the preferred cosmological frame - has been given comparatively little 
attention [1H7]. It has become accepted practice to simply remove the dipole from a sky map 
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before calculating the power spectrum |5]. This is based on the belief that the effect on Cg, the 
£th element of the angular power spectrum of the temperature anisotropy, (the central quantity 
in cosmological parameter estimation in the context of the canonical Lambda Cold Dark Matter 
model, aka ACDM) is proportional to /3 , where /3 ~ 10 -3 is our speed relative to the cosmological 
frame (as determined by the magnitude of that CMB dipole). In this context, the dipole is indeed 
the only multipole for which the Doppler shift due to our motion is significant. Nevertheless, it 
has been shown [9] that the second order (0(/3 2 )) Doppler effect noticeably alters the directions of 
the quadrupole multipole vectors, which characterize the shape of the quadrupole, despite the fact 
that it contributes negligibly to C2, the strength of the quadrupole. 

Meanwhile, several groups have independently examined (to order /3 2 ) the effect of Lorentz 
boosting the CMB, and determined that there are actually significant contributions to higher 
multipoles [IH7]. It has also separately been shown that simply masking a map induces correlations 
among nearby multipoles [10] . These two phenomena highlight a need for characterizing the effect 
of masking a Lorentz-boosted temperature field, since the process will undoubtedly mix higher 
multipoles of the true temperature field that have boost corrections. 

This paper is organized as follows: in Section III] we discuss the effect of boosting the multipole 



moments of the CMB, in Section III we derive an expression for the boosted pseudo power spectrum, 



in Section IV we show preliminary numerical estimates for corrections to the power spectrum and 



discuss their significance, as well as portions of the project that will be explored in future papers. 

II. THE EFFECTS OF BOOSTS ON THE MULTIPOLE MOMENTS OF THE 
ANGULAR POWER SPECTRUM OF THE CMB 

If an observer in the rest frame of the CMB (denoted S) measures a photon of frequency v 

arriving along a line of sight n, then an observer in another frame S' that is moving with respect 

to the CMB at velocity vv will measure the incoming photon to be arriving along a different 

line-of-sight n' with a different frequency v' . (Note that we will not concern ourselves here with 

any ambiguities in determining S associated with the existence of inhomogeneities, in particular a 

cosmological dipole, on the assumption that that dipole is 0(1O -5 ) smaller than the effects we will 

uncover.) The motion of the observer in S' thus induces two effects: a Doppler shift in the photon 

frequency and an aberration - a shift in the direction from which the photon arrives. These two 

effects can be seen explicitly in the relation between n and rV 

/ cos# + /3 \ „ n — cos#v 
11 = Vl + /3cos#y V+ 7(l + /3cos#)' ^' 



3 



where P = v/c and cos 8 = h ■ v. The change in observed frequency in S' is given by a simple 
Lorentz transformation 



v' = 71/(1 + /3cos<9) , 



(2) 



where 7 = (1 — P 2 ) 1//2 is the standard Lorentz factor. This angle 9 is related to the angle 9' , 
measured in the frame S' , via 



cos 0" 



1 cos # + /3 



(3) 



1 + P cos 6> 

Because of these two effects, we would like to boost the measured intensity, I(u, n) (the incident 
CMB power per unit area per unit frequency, per solid angle), and then relate the spherical 
harmonic coefficients of the intensity to the traditional temperature fluctuation coefficients. First, 
we will need the relation between the observed intensity in each frame 



I'(u, n') = 7 3 (1 + p cos 9f I(u, n). 



(4) 



Expanding both sides in terms of spin-weighted spherical harmonics ( s Y^ m ) and using the fact 
that dti = 7~ 2 (1 + P cos 6) dh |12j . we get the following expression for the boosted multipole 
moments, a'« m (y'\ in terms of the rest-frame multipole moments at m {v): 



t imi (v')= J dfl 7(l + /3cos6»)a fem2 (i/) s y£ 2m2 (n) s y/ imi (n). 



(5) 



We use the spin- weighted spherical harmonics to keep the expression completely general, so that the 
effect of the boost on the measured CMB spectrum can be explored for temperature fluctuations 
as well as polarization. Note that we have also implicitly chosen a frame where /3v is in the z 
direction. 

Expression §5§ can be expanded as a series in P by means of Eqs. Q and ^. The expression 
to order p 2 can be found in the appendix. To order P we find (consistent with [1]): 



4m (v') 



Psm 



2-1/ 



dv> 



where 



ssfm — 



at m [y ) 

1 {I 2 - m 2 )(e 2 - s 2 ) 
£ 2 (2£+ 1)(2£- 1) ' 



a^-l)m(^) (6) 
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Here, the multipole moments are the frequency dependent intensity coefficients. The more familiar 
temperature coefficients can be deduced from the coefficients above through the Stefan-Boltzmann 
law, which for fluctuations and a series expansion to G(/3) reads: 

8T(h) 1 61(h) 



T 4 I 

where 



(7) 



J(n) = / I(v,n)du (8) 



o 



is the incident CMB power per unit area per solid angle and To and Io are sky averages (monopoles) . 
Therefore, in order to get the temperature multipole moments we integrate Eq.Q over all frequen- 
cies (making use of the fact that the a!p m (v') are Planck distributed) and divide the result by four 
times the average intensity per solid angle |15| . To first order in /3 we have for the temperature 
a'lm (s = 0): 

a £m = a lm + /3Ci m a(e+l)m + /^ m a (^-l)m (9) 

where the conversion factor from temperature to intensity was absorbed into the multipolar coef- 
ficients and 

^ m = -(£-2) ^ +1)m , & m = (£ + 3)ottm- (10) 

We would now like to estimate the bias induced by the boost on the temperature power spec- 
trum. Before getting into the details of the calculation, it is first necessary to define two quantities. 
Throughout this paper we will make the distinction between the theoretical power spectrum, given 
by 

(a£ m a}, m ,) = 8 (ji 8 m>m ' Q, (11) 
and the measured power spectrum (also referred to as the power spectrum estimator), given by 



Ci = — 
1 21 



To estimate this bias we need to go second order in the expansion ^ since the C^s are quadratic 
in a£ m (see the appendix for details). The key point here is that, on the assumption of statistical 
isotropy of the unboosted ae m , the smallest j3 correction to the boosted power spectrum is given 
by @]: 

<C;>«C,(l + 4/3 2 + 0(/3 3 )). (13) 
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Note that the main effect is a rescaling of the spectrum by an overall amplitude 1 + 4/3 2 . Moreover, 
since (3 ~ 1(T 3 , the boosted power spectrum is essentially unbiased. We will now show that this 
conclusion may not be true in the case where the boosted spectrum is reconstructed from a cut-sky. 

III. MASKING EFFECTS ON THE BOOSTED ANGULAR POWER SPECTRUM 

When analyzing CMB temperature maps, it is common practice to mask regions of the sky that 
are believed to be contaminated. The region masked most often is the galaxy, where the observed 
temperature signal is known not to come from the surface of last scattering. In this case, we have 

a new expression for the measured temperature fluctuations: 

t 

AT (n) = E aimW(n)Y £m (n), (14) 
e>o m=-l 

where W(n) is a window function described by the mask. One can then decompose the left hand 



side of Eq.(14) and solve for a^ m in terms of a^ m . We then get the linear relation 

Q-£imi ^ ^2"i2 ^imifjm2 • (1^) 

Here Ki imi i 2m2 is a general kernel that contains all of the information about the window function, 
, w i imi Y^ imi (i\) . It is given explicitly by 

■(2*i + l)(2£ 2 + l)(2£ 3 + l)' 1 



4tt 




(16) 



h h is 
mi —rri2 ^3 

where the 3x2 matrices are the Wigner 3-j symbols. This can be used to find the pseudo power 
spectrum estimator: 



1 1 



Ci = — 

m=—t 

These Cgs are not unbiased estimators for the true C^s, but assuming statistical isotropy their 
expectation values are related by a real and symmetric mode coupling matrix, 

(4) = Y M ^(Ci 2 ) = Y, M hhCh, (18) 

where 

47r if 
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and We = l/(2£ + 1) ^ m |u;^ m | 2 is the power spectrum of the window function [TO] . 

We can therefore obtain an estimator for Ci by approximately re-writing the relation with {Ci 



replaced by its observed value on the sky Cf^ ~ Yl,i 2 We approximately invert Eq.(18) 

to obtain the true C^s from the pseudo ones. However, this prescription is complicated by Lorentz 
boosts. Since the masking procedure is unavoidably performed in the boosted frame, we no longer 



get a simple mode-coupling relation as in Eq.(18). Instead, we will now have a more general 
expression: 

1 h 
ft) = 21 + 1 ^ l a ^ m i' 2 

mi=— £\ 

= 21+1 ^ ^ ^ ^ a ^m2 a *hmz)^^imil2rn 2 ^} 1 m 1 tzrn- i ^ (^0) 
m\=—i\ l2m2 isms 

where the multipolar coefficients are given by Eq.([9]). Note that the kernels in this expression 
couple all the elements of the non-diagonal covariance matrix {a'i imi a} 2rn2 ) ■ As a consequence, 
linear j3 terms which were absent in Eq.([9]) will now contribute to the window function multipolar 
coefficients. 



Going back to expression Eq.(20), note that to first order in /3 there will only be coupling 
between £2 and £3 satisfying {£2 — 1| < £3 < £2 + 1- Therefore this expression can be re- written as 

1 h I2+1 

ft) = 2£ + 1 { a 'l2rn 2 a }[m 3 ) K hrn 1 l2rn 2 K } imi t 3 rn 3 - ( 21 ) 

m\=— 1\ f.2,rn,2 m 3 ,e 3 =l 2 — l 

When writing out each term in the sum over £3 and plugging in the expression for a!^ 2rri2 (see 
Eq.Q), we arrive at the following relation between the boosted pseudo-CVs and the expectation 
values of the CMB rest frame, C/s: 

ft) = (Mt&Ciz +Ne 1 g 2 Ce 2 + T^^C^+i + Q^C^-i) (22) 

where 

mi,m 2 

^i^a = 2£ Z 1 5^ Ke imi e2rn2Ke imi (e 2 +i) m2 £,t 2 m 2 > (^4) 

Qhh = 2£ +1 ^ ^' imifem2 ^imi(fe-l)m2^2m2 ' 
^ mi,m2 



Ignoring "fence-post" terms and re-indexing the sum, we can re- write Eq.(26) as 

ft) ~ ^ (AWa + A/^ 2 + P 4 ^ 2 -l + Qhi2+i) Ci 2 (26) 
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Clearly, the boosted pseudo power spectra in Eq.(26) are far from being unbiased estimators of 
the Ces. If we neglect the boost effects on the Cgs and only multiply by M.~ l to solve for the 
values of the C^s, then we will obtain an unbiased estimate since Af, V, and Q are /3-dependent. 



Furthermore, the bias in Eq.(26) produces more than just an overall multiplicative constant, since 
the matrices A4, J\f, V and Q are not in general peaked at £2 ~ l\. The correct reconstruction 
of the true C^s depends on the careful inversion of the sum of the four matrices above. Before 
carrying a complete numerical analysis, we must analyze the type of couplings that the matrices 
M, V and Q induce on the pseudo-power spectrum. This is the subject of the next section. 

A. Example: Equatorial Strip 

Before we calculate these couplings explicitly, we should mention two important points. First, in 
this section we will examine only the special case where both the boost and the normal to the mask 
are in the z-direction. This is of course unlikely to be correct. It is also likely to underestimate the 
magnitude of the effects we are examining. A thorough analysis that also includes the angle between 
the boost and the galactic plane as a free parameter is being carried out and will be reported on in 
a future work. Our point here is to note that (3 corrections exist regardless of this direction, and 
might in turn affect the reconstruction of the true Cgs if the boost is not correctly accounted for. 
Second, given that linear corrections will induce couplings of the form {£(1,2) > ^(2,1) ± 1) an d because 



of the restriction that £\ + £2 + £3 must be even in Eq.(16), the effect uncovered is a mixed-parity 
one. This means that f3 corrections are identically zero for masks which are even or odd functions. 
Real CMB masks do not possess well-defined parity, a fact which further points to the importance 
of the analysis being carried out here. 



To illustrate the type of coupling induced by the matrices ( 23p5 ), we have carried out a nu- 



merical analysis using a north-south asymmetric equatorial strip as a mask with different widths 
and degrees of asymmetry. Figure [T] show some plots of these matrices for £\ = 5 as a function of 
£2- Noticeably, the main difference between (23 25) and the matrix (19) is that the former assume 



both positive and negative values, whereas the latter is a positive-definite matrix. Regarding their 
amplitude, we have checked that matrices M, V and Q can be as large as 15% of matrix A4 for the 
range of £2 shown in Fig. [T] Including the f3 factor, this should amount to a correction of order 



10 -4 to Eq.(18). This may look smaller than the known 10 signal present in the off-diagonal 



terms of full-sky correlation matrices [HE]. However, we emphasize that here we are proposing a 
correction to the reconstructed true C^s, and not cross-terms of the correlation matrix as discussed 
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b 0.0005 - 



FIG. 1: The top figure shows the effect of moving the mask north by 5 degrees with respect to the galactic 
equator. Note that there are significant changes in contribution to the correction matrix elements as the 
mask is moved to 76-105 degrees (blue dots), 86-115 degrees (green triangles), and 96-125 degrees (red 
circles). The bottom figure shows the effect of varying mask size while keeping it only 1 degree north-south 
asymmetric for 76-105 (green triangles), 79-100 (blue dots), and 84-95 (red circles) degree masks. Here we 
notice slight changes in the matrix elements when the masks decrease in size, however the effect is very small 
when compared to the effect of moving the mask off-center. All of these values were generated with l\ = 5. 



in [7]; the former are not only less noise contaminated but also have a smaller cosmic variance. 
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IV. DISCUSSION AND FUTURE WORK 

Determining our direction and velocity with respect to the CMB rest frame is a fundamental 
quest for the standard cosmological model. While the main observable effects, noticeably the 
dipole, was already detected in the late 70's [13], higher multipolar distortions in the temperature 
field may yet be uncovered by future high-resolution data from the Planck satellite. In this paper 
we have shown that the CMB power spectrum may be systematically contaminated by this effect 
already at first order in the boost parameter if the boost effect is not taken into account when 
reconstructing full-sky data from cut-skies. We estimated this effect as a 10 -4 signal in the pseudo 
Cgs. Although this may seem smaller than the signal in off-diagonal full-sky CMB observables, 
only a careful reconstruction of full-sky data will be able to set this issue. We have also shown, 
using a simple strip as a mask, that the boost couples different multipoles in a way which is in 
sharp contrast to the couplings induced by the mask alone. Since the masking procedure is known 
to induce large-angle correlations, reconstructing the Qs correctly may have a sizable impact in 
tracing systematics and/or accounting for large-angle anomalies. 

In a companion paper we will carry a more complete analysis of these effects, including the angle 
between the boost and the galactic plane as a free parameter. This angle dependence will lead to 
further couplings between the boost and the mask which might be used as a further estimator of 
the boost direction. 

Since not all of the available information is contained in the Ci 's (due to the coupling between 
nearby modes of the a'g m ), we will define off-diagonal estimators of the covariance matrix. Addi- 
tionally, we will determine the likelihood of finding the direction of /3 from our boosted covariance 
matrix. Both of these issues have been discussed in part by [7j . We will also be looking into the 
effect of cut skies on reconstructing polarization estimators, which can be carried out starting from 
expression [6] and setting s = 2. 

Appendix A: Second order results 
1. Multipolar coefficients 

For completeness, we present here the expansion in the brightness multipolar coefficients to 
second order in the boost parameter. Similar expansions can also be found in [H [7] We now want 
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to expand Eq. [HJto 0(P 2 ). We begin by re-writing 



7(1 + cos 9) ~ 1 + cos 9 + -P 2 



and ae 2m2 (v) using 



/ «2 \ d 

at m {y) ~ 0^(1/) + ( /3 2 cos 2 - y - /3cos# j i/ —a tm {i/) + (p 2 cos 2 6) 

We will also expand s ^* m (n') about n to get 



,'2 



rf 2 



2 dz/ 2 



(Al) 



a lm {y') (A2) 



s Yg m (n')= s r/ m (n)-/3sin0(l + /3cosi 



dsYf m (n) , /3 2 sinW s Y/ m (A 



dcosO 



d cos # 2 



(A3) 



There is a general expression for writing the derivatives of a spherical harmonic in terms of other 
spherical harmonics (valid for £ > 1)[14|: 



sin ( 



d s Y em n (£+l)2_ m 2 



dcos# V (2£+ l)(2^ + 3) 



3 ^, m -(^+l)< 



£ 2 



(2£+ 1)(2£- 1) 



(A4) 



We also see that there will be some P cos 9 cross terms to deal with, which can be done using the 
formula 



cos 6 S Y{ 



s * £m 



(£+ l) 2 -m 2 



sYt+\,m 



sm 



s^f.m 4" 



sYl—l^ m 



(A5) 



' (2£ + 1)(2£ + 3) ^ * £(£ + 1) y (2£+l)(2£-l) 

Putting all this together and working to order p 2 , we arrive at an expression for the multipole 
moments measured in the moving frame as a function of those in the CMB rest frame, noting that 



the usage of A4 and A5 limits this expression to I > 2: 



Psm 



1 



l)m 



£(£ + 1) 

+ \s$m + 1)(* " 3) " 2(1 + 1)V'^ + 



/2 



d 2 



+ 



s 2 m 2 



, d v 12 d 2 \ 



e 2 {£ + 1) 2 V 2 " &v> 1 2 dz/ 2 y 



dz/ 2 j 2 



+ 4 ( 1 - i' 



dz/ 



ID l\ I / d 



2(^-i)-(^-2)i/— 



yd _ /2 d 2 



«(^+l)m( ly/ ) 



— >5.e 



Psm 



dz/ (£ +!)(£-!) 



l (£+2)m(V) 



+ \p 2 s^ + 2)mst ( l + l)m U* " 1) + 2^'~ + z/ 2 ^ 
+ ^/3 2 s 6ms^l )m (V + 1)(* + 2) - 2(* + + ^^2) 0(/-2)m(^) 



(A6) 
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where 

I { j2 _ m 2 ){i 2 _ S 2^ 

s ^ m ~ y £2(2* + l)(2*-l) ■ 1 j 

Note that this result differs slightly from the one shown in We have checked that they are the 
same up to some re-arranging. 



2. Expansion of the Pseudo Angular Power Spectrum 



Using expression (A6) and (20) we can show after a straightforward algebra that: 



£ 2 1 1% mim 2 

+ [/?& 2 m 2 (^2 + 3)1^4 mi^ 2 -lm 2 ] Q 2 -l}- (A8) 
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